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ABSTRACT 

- p  This  article  is  concerned  with  the  analysis  of  experiments  arranged  in  blocks  chosen  at  random. 

Estimates  of  the  polynomial  parameters  in  the  associated  response  surface  model  are  obtained  free  of 
blocks.  Tests  concerning  the  polynomial  and  random  block  effects  are  presented.  Furthermore,  the 
power  of  the  test  for  the  block  effect  is  obtained  using  a  certain  approximation  by  Hirotsu  (1979).  A 
numerical  example  is  given  to  illustrate  the  implementation  of  the  proposed  analysis. 
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1.  INTRODUCTION 

Model  fitting  in  response  surface  methodology  is  usually  based  on  the  assumption  that  the 
experimental  runs  are  carried  out  under  homogeneous  conditions.  This,  however,  may  be  difficult  to 
achieve  in  some  experiments.  For  example,  the  runs  may  be  obtained  over  a  period  of  time  or  in 
batches  among  which  the  experimental  conditions  can  vary  appreciably.  Such  an  extraneous  source  of 
variation  should,  therefore,  be  accounted  for  by  introducing  a  block  effect  into  the  response  surface 
model. 

Box  and  Hunter  (1957)  introduced  the  concept  of  orthogonal  blocking  so  that  estimates  of  the 
polynomial  parameters  in  the  mode!  can  be  obtained  free  of  the  block  effect.  This  concept  was 
developed  under  the  presumption  that  the  block  effect  was  fixed,  that  is,  represented  by  a  constant 
parameter  in  the  model.  In  this  case,  “the  effect  of  carrying  out  a  particular  trial  in  one  block  rather 
than  another  is  merely  to  change  the  expected  value  of  the  response  by  a  fixed  amount  which  depends 
only  on  the  particular  blocks  involved,”  as  was  stated  in  Box  and  Hunter  (1957,  p.  228).  Quite  often, 
however,  the  blocks  are  selected  at  random.  For  example,  the  blocks  can  be  batches  of  raw  material 
used  in  a  chemical  process.  In  this  case,  it  would  be  more  appropriate  to  regard  the  block  effect  as 
random. 

In  the  present  article,  the  analysis  of  a  blocked  experiment  will  be  discussed  under  the 
assumption  that  the  block  effect  is  random.  This  analysis  does  not  require  that  the  design  blocks 
orthogonally. 


2.  SOME  PRELIMINARIES 

Consider  a  model  of  order  d  (>1)  in  k  input  variables,  Xj,  y^,. . .,  x^.  The  experimental  runs 
are  arranged  in  b  blocks.  The  model  can  therefore  be  written  as 


y  —  /?gln  +  ?7  +  £» 


(2.1) 
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where  y  is  a  vector  of  n  observations  on  the  response,  /?q  and  0  =  (/3j,  0p )'  are  unknown 

parameters  associated  with  the  polynomial  portion  of  the  model,  7  =  (7^,  72,  ....  7^/,  where  7, 
denotes  the  effect  of  the  £—  block  (£  =  1,  2,...,  b),  and  r  is  a  random  error  vector.  The  matrices,  X 
and  Z  are  known  and  are  of  orders  nxp  and  nxb  and  ranks  p  and  b,  respectively.  The  matrix  Z  is 


actually  of  the  form 


?  -  ®  in,. 


where  n,  is  the  number  of  observation  in  the  £ —  block,  ln^  is  a  vector  of  ones  of  order  n,xl 
(£  =  1,  2,...,  b)  and  ©  denotes  the  direct  sum  of  matrices.  In  particular,  if  d  =  2,  that  is,  the  model  is 
of  the  second  order,  then  (2.1)  can  be  expressed  as 


yu  =  /?o  +  .£  Vui  +  £  Vui  +£.Vuixuj +  + 


,  u  =  1,  2,...,  n,  (2.3) 


where  denotes  the  u~  level  of  the  i~  input  variable  and  zu,  is  a  “dummy”  variable  taking  the 

value  one  if  the  u~  trial  is  carried  out  in  the  £~  block  and  zero  otherwise  (i  =  1,  2,...,  k; 
£  =  1,  2,...,  b;  u  =  1,  2,...,  n). 

Box  and  Hunter’s  (1957)  conditions  for  orthogonal  blocking  for  model  (2.3)  are 

£  xui  =  0,  i  =  1,  2,...,  k,  (2.4) 

u(£) 


ij  %  1,  2 . M*. 


y  x2.  =  ££  &  x2.  i  =  12  k 

Xui  n  ^ ,xui’  1  “  *'■■■'  K’ 

u(£)  u=l 


where  E  denotes  summation  over  the  £ —  block  (£  =  1,  2,...,  b).  Note  that  conditions  (2.4)  and  (2.5) 
u(£)  17 

presume  that  the  sum  of  the  elements  in  each  column  of  the  design  matrix  and  the  sum  of  the  cross 
products  of  the  elements  of  each  pair  of  columns  are  zero  (see  also  Khuri  and  Cornell  1987, 
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Section  4.7). 


If  the  block  effect  is  fixed,  that  is,  if  the  y^'s  in  (2.3)  are  unknown  constant  parameters,  then 
estimates  of  the  polynomial  parameters  in  the  model  can  be  obtained  in  the  usual  fashion  as  if  the  7^’s 

do  not  appear  in  the  model.  In  general,  for  a  model  of  order  d,  an  element  of  the  matrix  X  in  (2.1) 

6,  ^2  k 

can  be  written  as  xuj  x^  ...  xu£,  where  the  6j’s  are  nonnegative  integers  such  that  .  £^6j  =  S  and 

!<£<d.  In  order  for  a  design  to  block  orthogonally  we  must  have 


n  6j  &2 

£  xul  xu2 
u=l 


u£ 


lt)  =  0, 


It  follows  that 


«=1,  2,...,  b;  (2.7) 

6  =  1,2,...,  d. 


J 1  Jl  J k  _  n£  A  V61  *2  **k 

xul  xu2  **•  xuk  “  TT  2*  x ul  xu2  •••  xuk’ 
u=l 


which  can  be  written  as 


E 

u(£) 


£  =  1,  2,...,  b; 
6  =  1,  2,...,  d, 


*i 


E  3 

u(£) 


°1  °2 
ul  xu2 


xuk 


where 


£  =  1,  2,...,  b; 
S  -  1,  2 . d, 


(2-8) 


(2.8) 


[. 

is  a  design  moment  of  order  6 
moment  is  odd  if  at  least  one  6j, 


61  2  2 


kl  _  1  X^1 

J  —  n  2- xul 


u=  1 


ul  xu2  •  xuk 


(2.9) 


=  S  6;.  In  particular,  if  the  odd  design  moments  are  zero  (a  design 
i=l  1 

i  =  1,  2,...,  k,  is  an  odd  integer),  then  condition  (2.8)  requires  that 


. 

Condition  (2.8)  is  therefore  a  generalization  of  conditions  (2.4)  -  (2.6).  It  can  be  noted  that  the 
quantity  on  the  left  side  of  (2.8)  does  not  depend  on  subscript  £  and  is  thus  invariant  with  respect  to 
blocks  if  the  design  blocks  orthogonally.  In  this  case,  if  X  £  denotes  the  portion  of  the  matrix  X 
corresponding  to  block  £  (=  1,  2,...,  b),  then  from  (2.8),  the  lxp  vector, 
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{'  -  1^  ln*X*, 


(2.11) 


does  not  depend  on  subscript  l  and  is  therefore  the  same  for  ail  the  blocks. 

3.  THE  ANALYSIS  OF  EXPERIMENTS  WITH  RANDOM  BLOCKS 

Let  us  assume  that  the  block  effect  in  (2.1)  is  random  such  that  7  is  distributed  as  N^O,  <*7!^) 
independently  of  e,  which  has  the  normal  distribution  N^O,  <7^Inj.  Model  (2.1)  is  therefore  a  mixed 
model  since  0  is  a  fixed  parameter  vector.  In  this  case,  the  expected  value  and  variance-covariance 
matrix  of  y  are 


E(y)  =  /J0ln  +  (3-1) 

Var(y)  =  <r“Z?'  +  crfln  (3.2) 

=  er^-n'  (3-3) 

£=1  t 

where  Jn^  is  the  matrix  of  ones  of  order  n^xn^.  The  eigenvalues  of  Jn^  are  n*  and  0  of  multiplicity 

n.-l  (l  =  1,  2,...,  b).  By  the  Spectral  Decomposition  Theorem,  Jn  can  be  expressed  as 
*•  l 

Jn*  =  PfA*P'£,  i  =  1,  2 . b,  (3.4) 


where  A*  is  the  diagonal  matrix,  diag(n*,  0)  with  0  being  a  zero  matrix  of  order  (n*-l)x(n*-l),  and 
P*  is  an  orthogonal  matrix  of  the  form 


Note  that  for  l  =  1,  2,...,  b, 


l  =  I,  2,...,  b. 


(3.5) 


(3.6) 
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(3.7) 


9i9i  =  l«ri 
9l9't  - 

Formula  (3.3)  can  then  be  written  as 

Var(y)  =  <r|PAP'  +  4ln, 


where 


b 

P  =  ©  P. 
i=l  1 

b 

A  =  ©  A,. 

1= 1  1 

Now,  let  u  be  defined  as 

u  =  P'y- 

Then, 

E(u)  =  P'(/?01„  +  X0) 

Var(u)  =  +  <rfln. 

i=  1  1 


(3.8) 


(3.9) 


(3.10) 

(3.11) 


(3.12) 


We  note  that  the  elements  of  u  are  independent.  By  rearranging  these  elements,  it  is  possible  to 
partition  u  into  Uj  and  u0  of  orders  bxl  and  (n-b)xl,  respectively,  such  that 


»1  = 


b  -nr 

«=1  Ft  _ 


(3.13) 


u2  = 


b 


From  (3.1)  and  (3.3),  the  expected  value  and  variance-covariance  matrix  of  u0  are 


E(«2)  =  [  e%9't 

Var(u2)  =  o-?In_b. 


\X0 


(3.14) 


(3.15) 

(3.16) 
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3.1  The  Analysis  Concerning  the  Polynomial  Effects 


From  (3.15)  and  (3.16),  the  vector  u2  can  be  represented  by  the  model 


where 


u2  =  T-?  +  »?> 
1 =[«!,#■ 


(3.17) 


(3.18) 


and  »;  is  a  random  vector  distributed  as  N^O,  Since  X  is  of  full  column  rank  equal  to  p,  then 

T,  which  is  of  order  (n-b)xp,  is  also  of  rank  p  provided  that  n-b>p  and  the  matrix  S  =  [X:  Z]  is  of 
full  column  rank  (see  Appendix  A).  The  least  squares  estimator  of  0  is  therefore  given  by 


which  can  also  be  written  as 


0  =  (x'wx)  X'Wy, 


(3.19) 


where 

W  =  In  -  jJj  Jn£-  (3-20) 

Formula  (3.20)  is  true  because 

as  can  be  seen  from  (3.8).  In  particular,  if  the  design  blocks  orthogonally,  then  0  can  be  shown  to  be 
identical  to  the  least-squares  estimate  of  0  obtained  for  model  (2.1)  when  7  is  fixed.  The  proof  of  this 
assertion  is  given  in  Appendix  B. 

The  regression  and  residual  sum  of  squares  for  model  (3.17)  are 
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/  N-1 

u^t(t'T)  T,s2 

(3.21) 

.-1 

y'WXjx'WXj  x'Wy, 

(3.22) 

*2*2  -  SSReg 

(3.23) 

r  ,.,-11 

y  |_w  -  wx(x'wx)  x'vvjy. 

(3.24) 

It  is  interesting  to  note  that  SSp^  ls  in  fact  equal  to  R^/3  |  /?q,  7  j,  the  increase  in  the  regression  sum 
of  squares  due  to  adding  0  to  a  model  that  contains  the  intercept  0q  and  the  block  effect  y.  This  can 
be  clearly  seen  from  the  fact  that 

y'y  =  u'y 

=  UjUj  +  \J2U2 

b  r| 

~  «Sl  ^  +  SSReg  +  SSE’ 


where  r ^  is  the  total  for  block  l  (  =  1,  2,...,  b).  Hence,  the  total  sum  of  squares,  SS-pot,  for  model 
(2.1)  is 


SSTot  =H 


(3.25) 


b 

=  £  ... 
£=l  £ 


i 


+  SS 


Reg 


+  ssE. 


(3.26) 


But, 


(3.27) 
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is  the  increase  in  the  regression  sum  of  squares  due  to  adding  7  to  a  model  that  only  contains  the 
intercept.  This  is  the  usual  sum  of  squares  for  blocks,  SSgj^,  used  in  the  analysis  of  orthogonally 
blocked  experiments  with  fixed  blocks.  From  (3.26)  and  (3.27)  we  then  have  4 

SSTot  =  K?  I  *0)  +  SSReg  +  SSE*  (3'28> 

which  indicates  that  SSj^g  =  R^/?  |  0q<  7)-  In  particular,  if  the  design  blocks  orthogonally,  then 
SSpteg  =  |  @q)'  which  is  the  regression  sum  of  squares  for  the  polynomial  effects  obtained  by 

ignoring  the  block  effect  in  model  (2.1). 

When  replicate  observations  on  the  response  are  available  within  the  blocks,  the  residual  sum 
of  squares  can  be  partitioned  into  a  lack  of  fit  sum  of  squares,  SSpQp,  and  a  pure  error  sum  of  squares, 

SSpE.  The  latter  is  obtained  by  pooling  the  pure  error  sums  of  squares  from  the  blocks.  Tests 
concerning  the  polynomial  effects  can  then  proceed  using  SSpE  as  the  error  term  in  the  denominators 
of  the  F  test  statistics. 

3.2  The  Analysis  Concerning  the  Block  Effect 

The  sum  of  squares  for  blocks  given  in  (3.27),  namely  R^y  |  /?q),  is  not  appropriate  to  test  the 
hypothesis 

HQ:  <r^  =  0  (3.29) 

versus  the  alternative  hypothesis 

Ha*  9^  0, 

if  the  design  does  not  block  orthogonally.  This  is  because  the  expected  value  of  R^y  |  0^j  is  not  free 
of  the  fixed  polynomial  effects.  Since  model  (2.1)  is  basically  a  mixed  model,  a  more  appropriate  sum 
of  squares  to  use  to  test  the  block  effect  is  R^y  [  0q,  0*J,  which  represents  the  increase  in  the  regression 
sum  of  squares  due  to  adding  7  to  a  model  that  contains  the  intercept  0q  and  the  vector  of  polynomial 
parameters,  0.  The  expected  value  of  the  latter  sum  of  squares  is  free  of  all  fixed  polynomial  effects 
(see,  for  example,  Searle  1971,  Chapter  10).  More  specifically, 
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e[r(?  |  /fy  P)\  =  °7tr{?'  [in  -  X(X'X)  Vjzj  +  (b  -  1)4,  (3.30) 

where  X  =  [ln:  X].  Formula  (3.30)  follows  from  applying  formula  (79)  in  Searle  (1971,  p.  445). 
Furthermore,  under  Hq,  R.^7  |  /?q,  o\  has  the  chi-square  distribution  with  b-1  degrees  of  freedom. 
Thus,  an  appropriate  test  statistic  for  the  hypothesis  Hq  is 

r  "(?  I  »«■  <?)  ,3,,, 

F~(b-l)MSpr’  (3‘3  ' 


which  under  Hq  has  the  F  distribution  with  b-1  and  v  degrees  of  freedom,  where  u  and  MSpg  are  the 
number  of  degrees  of  freedom  and  mean  square  for  the  pure  error,  respectively. 

3.2.1  The  Power  of  the  Test 

Under  the  alternative  hypothesis  Ha,  R^7  |  /?q,  ^  j  <t\  does  not  have  the  chi-square 
distribution.  Instead,  R^7  |  $ J  is  distributed  as  I^AjVj,  where  the  V^’s  are  independent  chi- 

square  variates  with  one  degree  of  freedom  each,  and  the  A^’s  are  the  eigenvalues  of  the  matrix 

Q  =  (4 ??'  +  4ln)  S(S'S)  V  -  X(X'X)  V  ,  (3.32) 


where  S  =  [X:  Z].  See  Johnson  and  Kotz  (1970,  p.  151).  This  fact  can  be  used  to  obtain  approximate 
power  values  for  the  test  by  using  Hirotsu’s  (1979,  pp.  578-579)  approximation  (see  also  Khuri  1987. 
pp.  308-310).  More  specifically,  the  power  can  be  computed  using  the  approximation 


p(r2F».b-l.J“.)“p(rf.^h)  + 

[a/{3((+  2)(f+4)B(£,  £)}](! +  y») 


-(f+*)/2 


W  ^,f+2)(f+4)  +  24±i|±J 


2(f  +  t>)(f  +  4)  ,  (f  +  v  +  2)(f  +  v) 

i  .  //r  i.  \  i  O 


{l  +  i//(fh)}‘ 
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where 


F^  i  u  denotes  the  F  distribution  with  b-1  and  v  degrees  of  freedom  and  Fq  v  is  its  upper 
al00%  point,  B^,  denotes  the  beta  function  with  parameters  f/2  and  *//2, 

h  =  <r?(b  - 

f=[tr(fi)]2/tr(Q2), 

tr(Q)tr(Q3)  _ 

"  [tr(02)]2 

Formula  (3.33)  was  derived  from  formula  (3.5)  in  Khuri  (1987,  p.  309). 

4.  A  NUMERICAL  EXAMPLE 

The  yield,  y,  of  a  chemical  process  was  measured  at  various  levels  of  temperature,  catalyst 
concentration,  and  reaction  time.  The  coded  levels  of  these  input  variables  are  denoted  by  Xj,  x9,  and 
Xj,  respectively.  The  mean  response  was  represented  by  a  second-order  model  of  the  form 


Three  batches  of  raw  material  were  randomly  selected  and  used  in  this  experiment.  Each  batch  was 
only  large  enough  to  permit  a  maximum  of  eight  runs  to  be  made.  Some  variation  was  suspected  to 
exist  among  the  batches  since  they  were  received  at  different  times.  The  batches  were  therefore 
considered  as  blocks.  The  design  used  was  a  central  composite  design  consisting  of  the  three  blocks 
shown  in  Table  1.  This  design  is  rotatable  since  a,  the  value  of  the  axial  distance,  is  equal  to  1.682, 
which  is  the  fourth  root  of  8,  the  number  of  factorial  points  (see,  for  example,  Khuri  and  Cornell  1987, 
p.  118).  The  design,  however,  does  not  block  orthogonally  since  an  a  value  of  1.512  would  be  needed 
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for  this  purpose  as  can  be  verified  by  invoking  formula  (2.6).  The  observed  response  values  are  also 
given  in  Table  1. 


Using  formula  (3.19),  the  least-squares  estimate  of  /?,  the  vector  of  nine  coefficients  in  formula 
(4.1),  is  given  by 

0  =  (1.496,  1.145,  .382,  -2.734,  -1.701,  -1.037,  1.860,  -.958,  1.409)'. 

The  regression  and  residual  sum  of  squares  defined  in  (3.22)  and  (3.24),  respectively,  have  the  values, 
SSjieg  =  240.695,  SSg  =  7.4792  with  9  and  10  degrees  of  freedom,  respectively.  The  pooled  pure  error 
sum  of  squares  from  the  replicates  at  the  center  point  in  each  of  the  three  blocks  is  equal  to 
SSpE  =  2.0127  with  5  degrees  of  freedom.  Hence,  the  lack  of  fit  sum  of  squares  is  SSgQp  =  5.4665 
with  5  degrees  of  freedom.  Consequently,  the  lack  of  fit  test  statistic  has  the  value 


F  5.4665/5 
2.0127/5 


2.72, 


with  a  corresponding  p-value  of  .148  (see  also  Table  2). 


The  linear  and  quadratic  effects  associated  with  the  three  input  variables  are  displayed  in 
Table  2.  Their  sums  of  squares  form  a  partitioning  of  R^/?  j  0q,  j'j,  which  is  equal  to 
SSfleg  =  240.695.  On  the  other  hand,  the  sum  of  squares  for  blocks  adjusted  for  the  polynomial 
effects,  that  is,  R^7  |  0q,  0^  is  equal  to  121.4246.  Hence,  the  test  statistic  given  in  (3.31)  for  the 
hypothesis  Hq  in  (3.29)  concerning  the  block  effect  has  the  value 


F  -  121-4246 
2(2.0127/5) 


150.82, 


which  is  highly  significant  (the  p-value  is  .000034). 


5.  CONCLUDING  REMARKS 

The  conditions  for  orthogonal  blocking  impose  certain  constraints  on  the  settings  of  the  input 
variables.  Some  of  these  settings  may  not  be  feasible  from  the  experimental  point  of  view.  The 
analysis  described  in  this  article  can  be  applied  whether  the  design  blocks  orthogonally  or  not.  It  can 


be  conveniently  carried  out  by  using  standard  statistical  computer  packages  such  as  SAS.  This 
provides  the  researcher  with  more  flexibility  in  the  choice  of  design. 
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APPENDIX  A:  THE  RANK  OF  THE  MATRIX 
T  DEFINED  IN  FORMULA  (3.18) 


Let  r(-)  denote  the  rank  of  a  matrix.  Then 
r(T)  =  r(T'T) 


Now,  suppose  that  "  -ng/nf)J^  *8  not  c°lumn  rank-  Then  there  exists  a^O  such 

that  £  ^(in^  ~  -n^/n^J  Xa  =  0.  It  follows  that 


*»  =  [J,  (*»*/“«)]**- 


(A.l) 


[  f 

Xj:  X^:  ...  :  XM  ,  where  Xg  corresponds  to  block  l 
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(=  1,  2,...,  b).  From  (A.l)  we  then  get 


x.  =  [J,  Injc 


=  ZC,  by  (2.2), 


where  C  =  ^b)  an<*  ** t  ‘s  defined  by 


From  (A.2)  we  conclude  that 


C £  —  l  —  l,  2,...,  b. 


S[a':  -  C']  =  Q, 


where  S  =  [X:  Zj.  This  means  that  the  columns  of  S  are  linearly  dependent,  which  contradicts  the 
assumption  that  it  is  of  full  column  rank.  Hence,  the  columns  of  T  must  be  linearly  independent  and 
its  rank  equal  to  p,  the  rank  of  X. 


APPENDIX  B:  FORMULA  (3.19)  UNDER  ORTHOGONAL  BLOCKING 
Consider  formulas  (3.19)  and  (3.20).  If  the  design  blocks  orthogonally,  then  by  (2.11)  0  can 


be  written  as 


^(x'x-.fjf'^'y-Cj.yJe] 


-1  n«'(x'x)  /  „ 

=  (X'X)  Ip +  XV-^ 


-1  n(i'(x'x)  X'y  -  (  E  y  )( 

=  (X'X)  X'y  +  AK~-) 
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where  d  =  n^(x'x)  £. 


Now,  let  0t  denote  the  least-squares  estimate  of  0  when  y  is  ignored  in  model  (2.1).  This  is 
the  usual  least-squares  estimate  obtained  under  orthogonal  blocking  when  the  block  effect  is  fixed.  In 
this  case  0t  can  be  determined  from  the  equation 


,-lr 


r .  i 
00 

n  l'a 

& 

in 

h 

X'ln  X'X 

x' 

Hence, 


-1 


£r  =  (x'x)  X'y  + 


1  -  d 


which  is  identical  to  0  in  (B.l). 
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Table  1.  The  Design  Settings  and  Response  Values  for  the  Numerical  Example 


Uncoded  Variables 

Coded  Variables 

Response 

y  (oz) 

Temperature 

CC) 

Concentration 

(%) 

Time 

(h) 

X1 

x2 

x3 

170 

20 

3 

1 

-1 

-1 

75.505 

150 

26 

3 

-1 

1 

-1 

72.623 

150 

20 

7 

-1 

-1 

1 

69.377 

170 

26 

7 

1 

1 

1 

70.024 

160 

23 

5 

0 

0 

0 

70.171 

160 

23 

5 

0 

0 

0 

69.574 

160 

23 

5 

0 

0 

0 

70.179 

150 

20 

3 

-1 

-1 

-1 

57.202 

170 

26 

3 

1 

1 

-1 

68.801 

170 

20 

7 

1 

-1 

1 

69.686 

150 

26 

7 

-1 

1 

1 

69.460 

160 

23 

5 

0 

0 

0 

63.850 

160 

23 

5 

0 

0 

0 

63.787 

160 

23 

5 

0 

0 

0 

65.136 

176.82 

23 

5 

1.682 

0 

0 

75.300 

143.18 

23 

5 

-1.682 

0 

0 

72.281 

160 

28.05 

5 

0 

1.682 

0 

67.752 

160 

17.95 

5 

0 

-1.682 

0 

63.885 

160 

23 

8.36 

0 

0 

1.682 

72.752 

160 

23 

1.64 

0 

0 

-1.682 

72.278 

160 

23 

5 

0 

0 

0 

67.146 

160 

23 

5 

0 

0 

0 

68.253 
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Table  2.  Analysis  of  Variance  Table  for  the  Numerical  Example 


Source 

Degrees  of 

Freedom 

Sum  of  Squares 

(Type  I) 

F  Value 

p-Value 

Blocks 

2 

126.592 

Linear  regression 

X1 

1 

30.566 

75.93 

.00033 

x2 

1 

17.915 

44.50 

.00114 

x3 

1 

1.900 

4.94 

.07688 

Quadratic  regression 

X 
»— * 

X 

to 

1 

59.787 

148.52 

.00007 

xlx3 

1 

23.147 

57.50 

.00063 

X2X3 

1 

8.603 

21.37 

.00572 

0 

xi 

1 

52.429 

140.46 

.00008 

9 

x2 

1 

15.798 

33.95 

.00210 

x2 

x3 

1 

30.550 

83.80 

.00026 

Residual 

Lack  of  fit 

5 

5.4665 

2.72 

.148 

Pure  error 

5 

2.0127 
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